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Abstract 

Self-consistent factorization of two-body residual interaction is proposed for arbitrary density- 
and current-dependent energy functionals. Following this procedure, a separable RPA (SRPA) 
method is constructed. SRPA considerably simplifies the calculations and demonstrates quick 
convergence to exact results. The method is tested for SkI3 and SkM* forces. 
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I. INTRODUCTION 



Self-consistent nuclear models, as e.g. Skyrme-Hartree-Fock (SHF), are widely used for 
the description of nuclear ground state properties. They do also allow a description of 
excitation spectra. This is usually done within the random-phase approximation (RPA), 
which is meanwhile a textbook standard in nuclear physics [|l], 0, for a few recent applications 
see || |], [5] . Actually, these dynamical applications are up to now rather limited. The reason 
is that RPA implies the diagonalization of large matrices whose rank is determined by the 
size of the particle-hole lph space. The number of lph states grows huge for heavy and/or 
deformed nuclei. This limits the range of application for fully fledged RPA. The most 
applications are thus found for spherical nuclei. 

The RPA problem becomes much simpler if the residual two-body interaction is reduced 
to a separable form, i.e. 

v„ - vJi s ep) = 2 E ^z k z k ,, 

k,k'=l 

Z k = ^{p\Z k \h)a ] p a h (1) 

ph 

where {Z k } is a set of hermitian one-body operators. The factorization changes the rank of 
the RPA matrix from the number of lph configurations to the number of basis operators Z k 
and this reduces dramatically the computational effort. The success of a separable ansatz 
depends, of course, on a diligent choice of the operators Z k and the associated strength 
coefficients fi k y. 

Factorization of the residual interaction is widely used in nuclear theory but mainly 
within trivial schemes exploiting one separable term with an intuitive guess for the separa- 
ble one-body operator Z. The strength constant \i is usually fitted to reproduce available 
experimental data (see e.g. ||). Obviously, accuracy and predictive power of such sim- 
ple schemes are limited. Several improvements towards self-consistent schemes have been 



proposed during the last decades 0]-[12|. However, these schemes are not sufficiently gen- 
eral. Some of them are limited to analytic or simple numerical estimates others are 
not fully self-consistent in the sense that they start from phenomenological single-particle 
potentials flQ||-|ll|| or cover only particular effective forces ||12|| . 

In the present paper we propose a general self-consistent separable RPA (SRPA) approach 
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relevant to arbitrary density- and current-dependent functionals. The self-consistent scheme 
of Ref. H is generalized to the case of several separable operators. The operators are 
chosen to have maxima at different areas of the nucleus. This is crucial for an accurate 
reproduction of the residual interaction V Tes . A similar scheme has been successfully applied 
to the Kohn-Sham functional for description of collective oscillations of valence electrons 
in atomic clusters [I3|-[I5|. The nuclear case is more demanding since the SHF functional 
is much more involved. We will discuss the actual SRPA scheme for the case of SHF and 
present first successful tests for isoscalar E2 and isovector El giant resonances in 40 Ca and 
208 Pb using two typical Skyrme parameterizations. 



II. SRPA 



A. Separable operators 

In connection with nuclear density functional theory, it is preferable to sort the one-body 
operators Z k according to time parity. We thus decompose the set into hermitian time-even 
X k and time-odd Y k operators and the corresponding strength matrices K k k' and rj kk i. Then 
the separable Hamiltonian can be written as 

#SRPA = h + V^\ (2) 

TV sep 

K& P) = g ^(KwXkXv+muYkYv) (3) 
kk'=i 

where ho is the mean field Hamiltonian of the ground state, K kk > and rjkk' are symmetric 
matrices. Note that the operator product is considered to be separable at the Hartree level, 
i.e. for an arbitrary lph operator A we have 

[Z k Z k ,,A] = Z k ([Z k ,,A}} + Z k ,([Z k ,A]). (4) 

Furthermore, it is important to note for the following considerations that only commutators 
between time-even and time-odd operators have non-vanishing expectation values. That 
means that ([X k ,X«]) = ({Y k ,Yk>}) = while ([X k , Y k >]) ^ 0. 

The success of a separable ansatz depends very much on an appropriate choice of the 
expansion basis X k and Y k . Collective modes (vibrational, giant resonances) dominate the 
spectral strengths and it is intuitively clear that collective operators like the various multipole 
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moments should play the leading role in the expansion. We will denote them as hermitian 
time-even collective operators Qk (as e.g. a local multipole operator Qx^r) ~ r x (Yxfj, + Y^)) 
and time-odd conjugate momenta 

P k = i[H,Q k ] ph (5) 

where the index ph means the lph part of the operator. These operators can serve as 
generators of a collective motion by virtue of the scaling transformation 

|tf > = e d \), G — Y.imQh - iQkPk) (6) 
k 

where |) = \^q) is the ground state wave function, G is the generator, composed of various 
coordinate-like Qk and momentum-like i\ operators, Pk and qk are the corresponding (time- 
dependent) amplitudes, see e.g. [fTS], [T7|, |TB[. The dynamical deformation (||) induces changes 
of the mean-field Hamiltonian h which can be expanded in the linear regime for the small 

G as 

h ~ ho + [h , G] p h + [Kes> G]ph 

' — r: — ' (7) 

= ho + [H, G] p h- 

The dynamics of small-amplitude vibrations is governed by [H, G] p h and so we optimize the 
separable interaction by requiring that 

[H,G] ph = [H^\G] ph (8) 

as close as possible. Since both Hamiltonians contain ho in the same way, the requirement 
(|) is essentially 

[V ieS} G] ph =[V^\G] ph . (9) 
Applying (|9|) to generators Qk and Pk separately, we have 
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= J2x mkl {[X kf ,P k )), (10a) 

ibfe' 

[tL,4w = [Kl7\4u 

= EW6.<?*1)- ( 10b ) 



This shows that an optimal choice has to fulfill the condition 



X k e {[V zes ,P k ] ph j , (lla) 

ne {[K es ,4u}- (lib) 

We are free to assume any linear recombination among the sets {X k , Y k } and thus we may 
choose 

X k = [V res , P k ] p h, (12a) 
Y k = [V res ,Qk) P h- (12b) 



Substitution of Vjis^ into Q12a]) and ( |12b| ) has to give the identities. This results in the 



expressions for the inverse strength constants K k ,\ and r} w \ ($2 k K k l K kk' = hk', Y.kVjVkk' = 
S kk '): 

k~\ = ([X k ,,P k ]) = ([[V Ie8 , AU, A]), (13a) 

%\ = ([%,Qk]) = ([[Vres, Qk]ph, 4]>- (13b) 

The actual calculation of the commutators can be done either in wavefunction representation 



|19| , |20"fl or in terms of explicit lph matrix elements. The latter reads 

([X k , A]) = 2j2(p\Xk\h)(h\P k >\p), (14a) 

([y fe ,4']> = 2^( P |y fc |/ i )(/i|4b)- (i4b) 

Alternatively, the strength constants can be determined from the requirement 

([A, [v^s, A]]) = ([A, [v r t ep) , A]]), (isa) 

([4,[K es ,4']]> = ([4,[K { C T P) ,4]]). (isb) 



Substituting (|j) and ( |12a|) into ( |15a|) yields 



(fax,,]) = J2(i p ^ x k\)^m,p k ']) (i6) 



kk' 



and finally Eq. ( |13a|) . A similar procedure applies to the double commutator ( |15b|) . For 



a physical interpretation of the requirements ( |15a| ) and ( |15b| ) mind that the similar double 
commutators but with the full Hamiltonian (instead of the residual interaction) correspond 
to m.3 an d mi sum rules, respectively, and so represent the spring and inertia parameters 
11| in the basis of collective generators Qk and Pk- The condition ([13]) means that collective 



moments 7773 and m\ are exactly reproduced. 

The basic result of this section are Eqs. fll2f ) and ([0|) for the separable operators and 
strength constants, respectively. It is worth noting that these results can be also obtained 
pTH using the concept of nuclear self-consistency between the single-particle potential and 



density More detailed formulation of the basic double commutators in terms of the SHF 
functional will be given in section |IIE| and appendix |B[ The links to the vibrating potential 
model (VPM) |], [7| are outlined in appendix |A]. The choice of the generators Qk and Pk is 



yet open. An intuitive and pragmatic guess will be presented and discussed in sections [II F 
and 0. 



B. Conservation laws 



The forms ( |12|) are able to reproduce given conservation laws. Assume a symmetry mode 
whose generator P sym obeys by definition [H, P sym ] = 0. We simply have to include this 
mode into the set of generators, i.e. P sym G {Pk} together with its complement Q sym given 
by [H, Qsym] = —iPsym- One can prove that the separable Hamiltonian (||) does also fulfill 

[£ S RPA,Pym] =0. (17) 

From (1171) we have 



[/ig, P S ym] [^res; Psym] X ; 



sym- 
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Then 

[-^SRPA , -Psym] = [^0 ) Aym] + X k K kk i( [Xy , P sym ] ) 

k,k> 

— — X syra + X k Kkk'Kk'] 



,-1 
'sym 

k' 



which obviously yields Eq. (p!7[). Thus all symmetry modes are recovered if properly included 
in the ansatz. 

It is worth noting that in the case of a symmetry mode the operator of the residual 
interaction is 

X aym = — [^0) -Psym]) (18) 

i.e. can be also presented as the commutator with the single-particle Hamiltonian. In this 
case, the residual interaction serves to restore the symmetries violated in the mean field. 



The similar forces were proposed in [22] to restore translational and rotational invariance. 



C. Response functions 

One way to solve the SRPA equations is to compute directly a desired strength function 
S D {u) = ^\(* N \D\* )\\6(w-u N )+6(w + w N )) 

N 

= --a(([A 1 £>))) (19) 

where D is the operator of interest, and ojn are the wave function and energy of iV-th 
RPA state, L is the Liouvillian for small amplitude excitations. The latter is a shorthand 
for 

CA = [H : A] ph (20) 

where A is any lph operator. Similarly, we introduce the unperturbed Liouvillian Cq and 
the Liouvillian for the residual interaction V res by 

£ -4 = [h , A], v ICS A = [v res , A) ph . 
7 



For a compact formulation, we collect the operators D, {X k) k = l,...,7V sep } and {Y k , k = 
l,...,iV SC p} into the set {Z k) k = 0, 2iV scp } with X = D. The matrices of the strength 
constants are joined into the matrix ji of the rank 2iV scp + 1: 

( \ 0^ 

At = k • (21) 
v 77 ) 

Then we define generalized response functions 

1 



Rkk'(u) — ([Zk, 



uj + iT — C 
1 



Zk']), 

■z k ,]). 



(22) 
(23) 



uj + iT — C 

Note that the free response function can be written easily in terms of single-particle states 
p and h as 



ph 

1 



(24) 

^ + 2r-£ p?i UJ + lT + E ph ' 

where e p h is the energy of the lph-configuration and Hfc is the time-parity of the related 
operator Z k . 

A relation between full and free responses is established through the following steps: 

1 



Rkk' 



([z k , 



uj + iT-C - Vt p) 



Z k \) 



(25) 



i 



n=0 



uj + iT — C 



y(sep) 



uj + iT — C 



Z v \) 



= Ril + E R tl([z^ Wt p \ ^ 2 ]]><l + - 

- r(°) 4_ P(°) ,. 7?(°) 4. 

fclfc2 



This can be recollected to the final result 



JM*) = (« (0, (-)T^i^ 



MA./ (26) 

Note that this operator relation deals with matrices of the rank 2(2iV sep + l) (the additional 
coefficient 2 arises due to isospin). The strength function is finally 



S d (uj) = — 9{i2ooH}. 

7T 



(27) 
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The calculation of the free response function requires some caution with respect to time 
rity, see 
the pieces 



parity, see Eq. (p3|). The denominator in R^, mixes time parities and one has to disentangle 



1 ^1 0) M + 7^ 0) M 



uj + iT — Cq 



,(0)/ . UJ + iT 



(u + ivy - ci 

where Tl+ and connect operators of the same and opposite parity, respectively. One 
thus has to use nf for ([X,TZ^X]), ([Y,TZ^Y]), ([X,K<®D]), ([D,H^D]) and for 
([X,TZ^Y]), ([D,TZ^Y]). The mutually other combinations vanish. 



D. SRPA as eigenvalue problem 

As often done in RPA, one wants to determine the detailed eigenfrequencies and eigen- 
states of the response problem. The excitation phonon operator C\j for the mode N is 
determined from the RPA equations 

[H, C N ) p h = lo^C n . (28) 

The standard RPA scheme proceeds to recast that equations into the form of a matrix 
equation by expanding C^ N in terms of lph operators 

ph 

The involved matrices can grow huge in view of the large lph spaces required for suf- 
ficient convergence of the result. The eigenvalue problem is dramatically simplified by 
inserting the separable Hamiltonian (0). For compact notation, we recombine again the 
sets {Xk,Yk} — > {Zk, k = l,...,2A sep } and introduce the corresponding super-matrix of 
strengths (/ch',?|u') — > Mfcfe'- The RPA equations (pgj) then become 

Vkk>— , (30a) 

kk> N ~ 6ph 

(-) v^, (h\Z k \p)C Nik > 

c N P h = -l^Vkw — — — (30b) 

M UJ N + £ P h 
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with 



c N , k = ([i fc ,ci]> (31) 



,(±) 



Eqs. ( |30a| )-( P0~5D reduce the impressive number of unknowns c Nph to the much smaller set 



of unknonws Cjv.fc- Such reduction (which is possible due to factorization of the residual 
interaction) yields an effective RPA matrix with much smaller rank. 
As a next step, it is then convenient to transform C^,k to unknowns 

C-N,k = 2 H-kk'C-N,k' (32) 
k> 

which are directly connected with the amplitudes p k (t) and quit) of the scaling transforma- 
tion (H). For harmonic oscillations 

Qk(t) = qkCOs(u}t), pk(t) = p k sin(ujt) (33) 



and Cjv.fc = Qk and pk for time-even and time-odd cases, respectively. The unknowns (j32|) 
regulate self-consistently contributions of the scaling operators Qk and Pk into N-th RPA 
state. Every RPA state has it own set C^ ik an d pk) amplitudes. 

Finally, after standard algebraic steps, we obtain the SRPA equation for C^ jk 

^d kk ,(uj N )C N>k > = (34) 
k' 

with 

. / , n r (p|Zfcl fe )( fe l^fc'|p) £ pfa -i .oc-x 

rffcfc'(^v) = 2 2 _ 2 ^fefc" (35) 

ph £ ph U N 

Eq. (Q) has non-trivial solutions if 

det {d kk >(uj N )} = (36) 



which yields the RPA spectrum u^. Eq. ( p4|) is basically a matrix of the rank 4iV se p (after 
duplication of the rank because of the isospin) which is much smaller than the necessary 
number of lph states. Typically, we can deal with iV se p = 1 — 4 (see section [TIT]) as opposed 
to the rank of a conventional RPA matrix which often amounts to about 10 3 — 10 6 . 
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The normalization condition for RPA states is 



((C N ,Cj,)) = l (37) 

which amounts to 

£(i c £j 2 -i4^i 2 ) (38) 

ph 

Ex p sr^ < p\Z k \h >< p\Z k ,\h > e ph uj N 
CjV.fcCjv.fc/ 2^ , 2 _ 2A2 

fcfc' ph v P ft i ; 

1 - 9 

= 7 / CN,kCN,k'T\ dkk'i^N) = 2. 

4 tfc" dUN 

The last line of ([38|) represents the normalization condition in terms of derivatives of the 
RPA matrix elements. 

The reduced probability of i?A-transition from the ground state to the RPA state uj^ can 
be written as 

gr -> u N ) a~~mi — 7 vi ^ 69 > 

where D kk i{u N ) is the algebraic supplement of the matrix element dfcjfc'(wjv), 

^fc(Wjv) = 2^ 3 ~2 ( 40 ) 

and (h\f\\p) is the lph matrix element of E'A-transition. 



E. Connection with density functional theory 

Self-consistent models are often formulated in terms of a local energy-density functional 
rather than through an explicit Hamiltonian (see, e.g., the nuclear Skyrme-Hartree-Fock 
(SHF) method [23 or the Kohn-Sham scheme in electron systems [[24], |25|). So we will work 



out the above SRPA scheme in terms of a given energy functional 

E = E({J a (r)}) 

where 

J a (r) = <*|J (r)|¥) 
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constitute a set of local densities and currents associated with the functional. A simple 
example is the situation in atoms and atomic clusters where we deal only with the local 
density p(r) = |<p n (r)| 2 of valence electrons. The nuclear case is much more involved. 
The typical set for Skyrme functional includes, for both protons and neutrons, the local 
density p, kinetic density r, current j, spin-orbit density J, and spin density a i.e. 



The mean field Hamiltonian corresponding to the given energy functional is then deter- 
mined as 




(41) 



see in Appendix [B] the detailed expressions. 




(42) 



In the linear regime of small-amplitude deformation generated by a lph operator G, see Eq. 
(|), the densities can be decomposed as 




G 



(43) 



where J a (r) is the static ground-state density and 



4 G M) 



(*(t)|J (r)|«r(t))-(* |J«(r)|*o> 




(44) 





(45) 



We now decompose the scaling state as in Eq. (H), which yields the analogous decompo 
sition for the variation of density 




(46) 



k 



where, similar to Eq. (f|4j), 




(47a) 
(47b) 
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Every step proceeds much similar as for the general G because we are in the linear regime. 
Finally, we can specify the separable operators (|12a|) and ( |12b|) as 



n (r) =E/^-w fJ . ( g, M J ? t W- (48b) 



a/3 

and the strength matrices ( |13a| ) and (|13b| ) as 

«S =X/*drVf> (r) JJ £±-^- ) J^ (r') , (49a, 

=E/^'^(r) , Ja( g g(r0 ^-(r-) • (49b) 

The latter equations demonstrate the symmetry of the strength matrices and r^/. Note 
the subtle difference between separable operators and strengths. The double integrals look 
much the same. However, in Eqs. ( }48a| ) and ( [48 b| ) , they contain the current operators J a (r) 
and yield one-body operators as a result of the integration. Instead, Eqs. (|49a|) and 



do not include operators and so yield strength coefficients as c-numbers. For more details 
see Appendix B. 

F. Choice of collective generators 

As was mentioned above, the proper choice of generating operators Q/c(r) is crucial to 
achieve good convergence of the separable expansion (|3]) with a minimal number of separable 
operators. The choice is inspired by both physical and computational arguments. It should 
be simple and universal in the sense that it can be applied equally well to all modes and 
excitation channels. The main idea is that the generating operators should explore different 
spatial regions of the nucleus, the surface as well as the interior. This suggests that the 
leading scaling operator should have the form of the applied external field in the long-wave 
approximation, for example, Qfc=i( r ) = rA (Xv(^) +h.c). Such a choice results in separable 



operators (48a) and ( |48b|) which are most sensitive to the surface of the system. The first 
excitation modes of nuclear collective motion tend to vibrate predominantly in the surface 
region. As a result, already with this one separable term we obtain a quite good description 
of giant resonances. The detailed distributions depends on a subtle interplay of surface and 
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volume vibrations. This can be resolved by taking into account an interior of the nucleus. 
For this aim, the radial parts with larger powers r A+p Yx^ and spherical Bessel functions can 
be used, much similar as in the local RPA [|1^, |20[ . This results in the shift of the maxima of 
the operators ( |48ap and Q48bp to the interior. Exploring different conceivable combinations, 



we have found a most efficient set of the generators as 

Q fc (r) = R k (r)(Y x ^n)+ h.c.) , P k = t[H,Q k ] (50) 

with 

\ * , (5l) 

Jx (q k x r), A; = 2,3,4 

<?a = a fc"5^> a 2 = 0.6 ,a 3 = 0.9 ,a 4 = 1.2 

where i?diff is the diffraction radius of the actual nucleus and z\ is the first root in j\(z\) = 0. 
The separable term with k — 1 is mainly localized at the nuclear surface while the next 
three terms are localized more and more in the interior. This simple set seems to be a best 
compromise for the description of nuclear giant resonances in light and heavy nuclei. 

One may argue that one needs to explore not only different spatial regions but also the 
different types of the operators, associated with the variety of densities and currents (e.g. 
kinetic, spin-orbit). We have checked this. And indeed, this allows a perfect fine tuning of 
the spectra. On the other hand, it adds more separable terms. The set (p0|) -(pl|) with few 
and purely local generators is in our opinion the best compromise between quality of the 
results and the expense of the calculations. 



III. RESULTS AND DISCUSSION 



We test SRPA for Skyrme forces by comparison with standard full RPA calculations. As 
the test cases we consider the isovector dipole and isoscalar quadrupole resonances in the two 
doubly magic nuclei 40 Ca and 208 Pb. Most tests are performed with widely used Skyrme force 
SkM* pOj. The convergence of the SRPA expansion is found to be only slightly dependent 
on the actual force, which is demonstrated for the case of the recent parameterization SkI3 
||28|| . The full RPA calculations are done in a large basis which covers the lph excitations 
from all hole states to all particle states up to 30 MeV above the Fermi energy. SRPA easily 

14 



allows to extend this basis but we restrict the lph summations in SRPA to precisely the 
same active space as used in full RPA to have one-to-one comparison. The comparison is 
done in terms of the strength function <Sd(u;) using a smoothing parameter of T = 1 MeV. 
Specifically, the Lorentz weight function was used. The value of the smoothing is quite 
realistic because experimental spectra are typically washed out by similar widths. 

Figure 1 compares SRPA and full RPA results for SkM* forces. Three stages of SRPA 
expansion with increasing number of separable terms are shown: only one term k = 1, 
two terms k = 1,4, and all four terms according to the labeling (f)T|) . Mind that we use 
in the second set the operator with = 1.2 since then the forces cover the extremes of 
both surface (k=l) and volume (k=4) sensitivity. The first stage with only k = l provides 
already very good description of the isoscalar E2 resonance in both nuclei. Even the low- 
energy quadrupole mode in 208 Pb is well reproduced. The influence of further terms is very 
small. The isovector El resonance seems to be more demanding. The first stage correctly 
describes the average position and width of the resonance but cannot fully cope with the 
detailed fragmentation pattern. Here one sees clearly the systematic improvement achieved 
by the further terms. The set with four terms provides a satisfying agreement. Further 
improvement is possible by adding more terms but the extra effort seems to be unnecessary. 
The agreement achieved is already more than sufficient for any practical calculations. One 
may even be content with a smaller sets for large scale exploratory studies. 

Although the smoothing with 1 MeV provides the typical realistic shapes of strength 
distributions, it is worth having a quick glance at spectra with higher resolution. This 
is done in Fig. 2 drawn with smoothing T =0.1 MeV. The figure provides the detailed 
comparison of SRPA and exact RPA results for the most complicated cases of the isovector 
El resonance. Though we have not a perfect coincidence (which may be partly excused by 
numerical noise in the full RPA calculations), the results are indeed extremely close. This 
confirms the reliability of SRPA even at this finer scale. 

Figure 3 demonstrates the quality of SRPA for another Skyrme force, namely the pa- 
rameterization SkI3. It shows the toughest test case, the isovector El resonance in 208 Pb. 
The full set k=l-4 is used. Again we see a satisfying reproduction of this "worst case". The 
other test cases (not shown here) perform even better. We also checked a few other forces 
and always found the same quality. 

In fact, the spectrum from SkI3 looks much more realistic than that obtained from SkM* 
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(see lower right panel in Fig. 1. The pronounced right shoulder in SkM* is, to a large 
extend, an artefact from the restricted lph space used in these calculations. SRPA offers 
the opportunity to remove that restriction and to perform calculations in the full space of 
excitations, limited only by numerical resolution of the underlying coordinate space grid 
(with particle energies at least up to 600 MeV). Figure 4 compares SRPA results in the 
restricted and unrestricted spaces. It is obvious that the consistent inclusion of lph states 
above 30 MeV serves to diminish the unnaturally high shoulder at 16 MeV. This is a clear 
hint that RPA calculations converge slowly with the size of phase space. The SRPA allows 
to go much further in that respect. 

Altogether, the tests show that the proposed set of four operators suffices for a satisfying 
reproduction of full RPA. The number "four" has to be taken with a grain of salt. We 
have four operators of the type Qt for protons and another four for neutrons. The set is 
complemented by the conjugate operators Pk = i[H,Qk\- This amounts to 16 operators in 
total and the computation of the strength function requires inversion of 17 x 17 matrices. 
The full RPA calculations, on the other hand, used matrices up to a rank of 400 and more. 

IV. CONCLUSIONS 

A novel self-consistent separable approach to the random phase approximation (SRPA) 
is proposed. SRPA employs a systematic expansion of the exact residual interaction into a 
sum of weighted separable terms (i.e. products of one-body operators evaluated at a Hartree 
level). For both, the weights and the separable operators, compact analytical expressions 
in terms of collective generators are derived in a self-consistent manner. The form of the 
separable residual interaction is optimized by the condition that the vibrating mean field as 
generated from collective flow reproduces the exact result. SRPA is formulated in a general 
way such that it can be applied to arbitrary energy functionals depending on various local 
densities and currents. As open point remains the choice of the relevant collective generators. 
Here we exploit experience gained in nuclear fluid dynamical models (local RPA) and use 
an efficient mix of surface and volume modes. 

The general SRPA scheme is specified for the particular case of the Skyrme energy func- 
tional. The convergence of the separable expansion and performance of SRPA for giant 
resonances in light and heavy nuclei is investigated. The results are very encouraging. 
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SRPA proves to be a very good approximation to exact RPA and it it requires much less 
numerical effort. In particular, it allows to work with huge particle-hole spaces, giving thus 
the chance for unrestricted RPA calculations, even in nuclei with broken symmetries 
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APPENDIX A: COMPARISON WITH VPM 



It is interesting to compare SRPA with the Vibrating Potential Model (VPM) |1|, 0] whose 
different versions are widely used in nuclear structure theory (see, for example, || |9], [It], [EL 



29[). VPM deals with the mean field 

h = T + U (Al) 

where Uq is purely local. Moreover, Uq is often a density-independent phenomenological 
potential (oscillator, Nilsson, Woods-Saxon, ...). The vibrations of the system with the 
mean field ( |A1| ) are characterized by the collective deformations q k and in the framework of 
VPM the vibrating potential is written as 

dU 

U = U + 5U, 5U = S2—5q k = -tJ2x k 5q k (A2) 
where X k —i^- are single-particle operators in the r-representation, participating in VPM 



separable residual interactions (see e.g. [p9|0. Usually the collective momenta P k are intro- 
duced 

P k = i^-. (A3) 
dq k 



This allows to write 



X k = i^ = [U,P k ] 
oq k 



X k = j drp(r)[U, P k ] « J drp(r) < [U, P k ] > (A4) 
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where X k is the operator corresponding to X k and acting in the Fock space. VPM operators 
(|A4|) should be compared with the SRPA operators ( |1 2a| ) or (|48aj) . For the first sight Eq. 
( |A4l ) contradicts to ( |12a| ). However, there are, in fact, more similarities than differences. 
The equivalence of (|A4|) and ( |48a| ) can be demonstrated in the case of the simple energy 
functional 



E = \t j dv p{vf 



This functional gives the following mean field potential 



U = t p. 



(A5) 



(A6) 



VPM expression ( |A4| ) gives 



X k = t drp(r) < [p{r),P k ] >, 



(A7) 



while SRPA expression (|48a|) provides 



Xh 



drdr' 'p(r) 



6 2 E 



5p(r)5p(r 
t /drp(r) < [p(r), A] > • 



-<[p(r'),P k }> 



(A8) 

Thus there is no difference between VPM and SRPA for simple two-body forces following 
from the functional (|A5|) . Differences come up with non-linear terms in the mean- field 
potential. From this point of view SRPA seems to be more general. 



APPENDIX B: DETAILS OF SHF 



A most widely used energy functional for nuclear structure calculations is the Skyrme 
functional. It was originally proposed in P^] and made its breakthrough with the first fine- 
tuning of [27]. Since then it has been steadily further developed and constitutes today a 
very reliable model for nuclear structure and excitations. The functional involves the set of 
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local densities and currents ( [41]) which read in general 



J<?( 



$^n( r Vn( r ), 
^V^,(r)-V^,(r), 

n£q 

-^<W(V x <3>„(r), 

n£q 

-^E[<W V ^W-V^(r)^(r)] 



The associated operators are 



n€q 



Pq( r 



o- g (r 



iV(r'-r), 

n^^f - r)p, 
fl g 5(r' - r)px<j, 

in,{p,^-r)} 

ng5(r' — r)a 



where fig is the isospin projector to q G {p, n}. 



We use here the Skyrme functional in the form [pO 

E = J dr(£ kin + £ Sk (pq,Tq,a q ,i q , Jg)) 



^kin 
Eq 



+Ec(p P )-E c 

h 2 



2m 

-,2 



-T. 



e 2 f 1 
— J dr dr p p (r)———p p (r) 



3 o / 3 



(Bl) 
(B2) 

(B3) 



<? <? 
-|pAp + P« A Po - &4 (pV • J + a • (V x j)) 

-K • J «) +v(vx j ff )i - w 2 -^E4 



(B4) 
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The sum Yl q runs over proton or neutron species. A density without index means total 
density, e.g. p = p p + p n and similarly for the others. The parameters hi and b[ used in the 
above definition are chosen to give a most compact formulation of the energy functional, 
the corresponding mean-field Hamiltonian and residual interaction . They are related to the 
standard Skyrme parameters £j and x» |27| by 



bn 


— tn(l + -Xn) 




b'o 


= t (| + x o)? 




h 


= \ + 


+ t 2 (l + \x 2 )\ 




= i [tiii + xx)- 


-t 2 (l + x 2 )] , 


h 


= | [3*1 (1 + \x x 


)-t 2 (l + \x 2 


b> 2 


= |[3ti(| + xi) 


+ t 2 (l + x 2 )] 


h 


= It 3 (l + Ix 3 ), 




b' 3 


= 4^3(2 + ^3), 




bi 


= b' A = |*4. 





(B5) 



Various versions exist for the spin-orbit term. The J 2 term is sometimes considered, some- 
times not. Thus one has the choice 

^ = -je{hxi+t 2 x 2 ), 

(B6) 

b' 4 = * 2 ) 

or 

b A = 0, 64 = 0. (B7) 

The pV- J term is always taken into account. Two options exist here. Conventional Skyrme 
forces use b' A = 64 while a recent extension allows for more flexible isospin dependence 
considering 64 and b' 4 as independent parameters E3]. The center-of-mass correction is 



-^cm = (-P C m)/2 m A often approximated by its diagonal terms E cm = (J2iPl) /2mA. The full 
correction is usually subtracted a posteriori and the approximate correction is expectation 
value of a one-body operator. In any case, there is no contribution to the residual interaction 
for RPA. 

First variation with respect to J a yields the mean-field Hamiltonian. It can be found 



in many publications, e.g. [20|. The step of interest is the second variation yielding the 



residual interaction. We obtain here 
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5 2 E 



^(r)^(r') " I 60 + 63 3 P " h "^^ P \ P «- b3 Y P ~ b2A ) 6{r ~ " } 

~(b' + b'^p* + b^-pf^p, + b' 2 A)6 qq ,6(r - r') 

+M P ,<r^-^ 2 (^P; 2/3 ^(r-r')), 



5 2 £ 



5p ? (r)5v(r' 
<5 2 £ 



<Jp,(r)<JJ,/(r' 
<5 2 £ 



5J i/g (r)5J i/ y(r / 
5 2 E 



5 2 E 



M(r-r , )-6i<J(r-r')V, 
-(fe4-^)v5(r-r'), 



&4 + &4 ( Vj ( K r - r ')<W> 

2& 1 -2&' 1 v)<H r - r ')<W, 

(& 4 + ^vjv^ - r ')w M 



<5j^(r)5fVg/(r') 

where z/, z/ e {x, 2} labels the spatial vector components and e w ij, is the antisymmetric 
tensor. All combinations of derivatives not listed above vanish. 

It is instructive to write the response currents in detail. The non-vanishing contributions 

are 

pf(i-) = < [p q ,P k ] > 



n£q 



7f (r) = <KPk}> 



^(v<-vA^-vA<-v^), 



n€<j 



Jf(r) = <[J ff ,P fc ]> 



n€<j 



j?*(r) = <[j 9 ,g fc ]> 



9 1 Y (yl^Qkifn - Q k ip* n Vcp r 



(7^(0 = < [& q ,Qk\ > 



n£q 
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FIGURE CAPTIONS 

Figure 1: Isovector El and isoscalar E2 giant resonances in 40 Ca and 208 Pb calculated 
with SkM* forces in full RPA and SRPA at three levels of expansion as indicated. The 
strengths are weighted by Lorentz function with the averaging parameter T = 1.0 MeV. 

Figure 2: Isovector El resonance in 40 Ca and 208 Pb calculated with SkM* forces in full 
RPA and SRPA with the complete set k — 1 — 4. The responses are weighted by Lorentz 
function with the small averaging parameter T = 0.1 MeV. 

Figure 3: Isovector El resonance in 208 Pb calculated with SkI3 forces in full RPA and 
SRPA with a set of four operators. The responses are weighted by Lorentz function with 
averaging parameter T = 1.0 MeV. 

Figure 4: Isovector El resonance in 208 Pb calculated in SRPA with the standard set 
of four operators and using the SkM* force. Two different phases spaces are compared: 
full line = unrestricted space of excitations, dashed line = space restricted to ph energies 
< 30 MeV as in previous figures. The responses are weighted by Lorentz function with 
averaging parameter T = 1.0 MeV. 
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